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Abstract 

The masses and decay constants of the hght vector mesons p/uj, <j) and K* are 
studied within a ladder-rainbow truncation of the coupled Dyson-Schwinger 
and Bethe-Salpeter equations of QCD with a model 2-point gluon function. 
The approach is consistent with quark and gluon confinement, reproduces the 
correct one-loop renormalization group behavior of QCD, generates dynam- 
ical chiral symmetry breaking, and preserves the relevant Ward identities. 
The one phenomenological parameter and two current quark masses are fixed 
by requiring that the calculated /t^, m-Tr and are correct. The resulting 
fx is within 3% of the experimental value. For the vector mesons, all eight 
transverse covariants are included and the dominant ones are identified; the 
complete angle dependence of the amplitudes is also retained. The calculated 
values for the masses rup, and tuk* are within 5%, while the decay con- 
stants fp, ftf, and fx* for electromagnetic and leptonic decays are within 10% 
of the experimental values. 
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I. INTRODUCTION 



A realistic description of vector mesons at the quark-gluon level is an important ele- 
ment in advancing our understanding of hadron dynamics and reaction processes at scales 
where QCD degrees of freedom are relevant. They are easily produced as decay products 
in electro-excitation of baryon resonances and also as precursors to di-lepton events in rel- 
ativistic heavy-ion collisions. Flavorless vector mesons couple directly to the photon and 
play an important role in the phenomenology of electromagnetic coupling to hadrons. This 
is exemplified by the general phenomenological success of the Vector Meson Dominance 
model which assumes that the electromagnetic current is saturated by the vector mesons. 
The ground state vector mesons, being spin modes, sample the qq bound state dynamics 
in a way that is complementary to that of the ground state pseudoscalar (PS) mesons. 
They also explore quark and gluon confinement since the vector masses are greater than 
the sum of typical constituent quark masses. Mesonic strong decays such as p tttt and 
(j) — > KK, radiative decays such as p ^ 7r7 and K*^ — > K'~^, and electromagnetic decays such 
as p° — * e^e~ and (j) — > e^e~ can probe aspects of the underlying quark-gluon dynamics that 
are complementary to what is learnt from PS mesons. 

The PS mesons, especially the pion and kaon, have for a long time been a major focus 
of attempts to understand the internal structure of hadrons from nonperturbative QCD. 
Chiral symmetry provides an assistance in the PS case that is not available to other mesons 
like the vectors we discuss in this paper. In the chiral limit of massless QCD, the dynamical 
breaking of chiral symmetry generates masses for the light quarks that are consistent with the 
empirical constituent masses deduced from the hadronic spectrum. The lowest PS mesons, 
which would be massless in this limit as dictated by the Goldstone theorem, acquire a mass 
through the explicit breaking of chiral symmetry induced by the current quark masses. 
This phenomenon dominates the systematics of the ground state PS octet and provides 
chiral Ward identities that relate dynamical quantities in a way that simplifies somewhat 
the task of modeling low energy QCD. For example, the axial Ward identity dictates that 
the chiral limit Bethe-Salpeter (BS) amplitude for a pseudoscalar qq bound state in the 
dominant 75 channel is given by Bq{p'^)/ fp where Bq is the scalar part of the chiral quark 
self-energy, and fp is the meson weak decay constant. Consideration of the symmetry 
breaking effect of current masses within the PS bound state dynamics leads to an exact 
formula for mp the PS meson mass One corollary is the Gell-Mann-Oakes-Renner 
relation at small current quark masses where mp (x y/rriq] a second corollary is the behavior 
mp (X mq for heavy quark PS mesons [0. With these aids, an efficient and qualitatively 
useful phenomenology for observables and other quantities associated with the pion and 
kaon can be produced without explicit solution of the bound state Bethe-Salpeter equation 
(BSE). The only dynamical input required is the dressed quark propagator as defined by 
the quark Dyson-Schwinger equation (DSE). 

The study of hadronic processes is often facilitated by parameterizing DSE quark propa- 
gator solutions into analytic forms with a few parameters readjusted to accommodate chiral 
observables 0. With such an approach, m^/i^ and /vr/i^, the charge radii t^jk and form 
factors Fj./KiQ'^), and the vr — vr scattering lengths have been well described This ap- 
proach has also been successful in studies of coupling constants and form factors for processes 
such as 7r° — 77 [§], and •yn ^ nn f^. The incorporation of vector mesons has been hin- 
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dered by the lack of a symmetry-based means of obtaining approximate dynamical insight 
without direct solution of the vector BSE. There is a conserved current and a correspond- 
ing vector Ward-Takahashi identity linking the longitudinal vector vertex with the quark 
propagator. However, this Ward identity does not constrain the transverse vector meson 
BS amplitude and therefore purely phenomenological vector BS amplitudes have often been 
used to study processes involving vector mesons in this framework. Successful applications 
of this type include the decays p —>■ tttt and p ^ tc'j [^,|^ and diffractive electroproduction of 
vector mesons 0. Uncertainties concerning contributions to hadronic observables from the 
neglected covariants in the PS meson BS amplitudes beyond the canonical 75 were not ad- 
dressed until recently Several studies also incorporated some aspects of vector 
BSE dynamics and were able to make a crude assessment of the role of sub-dominant covari- 
ants of the rho [|1^] in p ^ vrvr and p — 717 but at the expense of using a separable Ansatz 
for the BSE kernel. A recent study of heavy meson decays also employs phenomenological 
BS amplitudes for vector mesons ||^. 

The persistent outcome of the above studies is that soft observables associated with the 
pion and kaon are consistently and naturally described in terms of the momentum dependent 
quark self-energy from realistic solutions |]13[ of the quark DSE. Important to the success 
of this approach are the features of quark confinement and dynamical chiral symmetry 
breaking that are implemented through a strong enhancement in the infrared behavior of 
the effective quark-quark interaction (or effective gluon 2-point function |T^) in rainbow 
approximation. In the light pseudoscalar sector, the most comprehensive and quantitatively 
reliable study to date |]10| involves direct solution of the bound state BSE in conjunction 
with quark DSE solutions for propagators. That work represents the development of an 
appropriate phenomenological representation for the infrared structure of the gluon 2-point 
function in conjunction with a bare quark-gluon vertex so that the DSE solution for the quark 
propagator exhibits dynamical chiral symmetry breaking as well as confinement [l^ and, 
through the BSE, produces a good description soft pion and kaon observables. One of the 
aims of Ref. |]10[ was an exposition of the detailed numerical consequences of the constraint 



provided by the axial vector Ward-Takahashi identity (AV-WTI) upon PS meson dynamics. 
This constraint is formally assured by the coordinated rainbow-ladder truncation of the DSE- 
BSE complex of DSEs. The inclusion of all possible covariants for the BS amplitude was 
found necessary to numerically preserve the constraint and to obtain quantitatively accurate 
observables. Of general importance are two other advances represented by Ref. [0 . Firstly, 
since the ultraviolet structure of the model is equipped with the one-loop renormalization 
group properties from perturbative QCD, it is a realistic and covariant hadron model that can 
be unambiguously evolved in scale. Secondly, it is produced by well-defined truncations of 
the QCD equations of motion (DSEs), and thus can be systematically improved by including 
higher-order corrections to the quark- ant iquark scattering kernel. 

The extension of the DSE approach to vector mesons is explored here. Solution of 
the vector BSE is more difficult than in the PS case because of the significantly larger 
number of covariants that must be investigated and also because the higher masses produce 
a larger domain of the quark complex plane that must be sampled. This latter issue 
was avoided in a previous work that made an extensive study of the meson spectrum 
from the ladder-rainbow truncation of the DSE-BSE system. In that approach, a derivative 
expansion of the quark self-energy was used to infer the behavior away from the real axis 
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and some attempt was made to estimate the resulting error. The imphcations for quark 
confinement in that approach are unclear. One of our aims here is to generate vector meson 
BS amplitudes without compromising the analytic structure in a way that may impair the 
subsequent explorations of meson decays and form factors. These amplitudes can then be 
used to calibrate and guide approximate representations that simplify the study of hadronic 
interactions. 

In this paper we calculate the ground state vector mesons p/cu, K* and cf) in the DSE- 
BSE approach, using the ladder-rainbow truncation. The effective quark-quark interaction 
is fixed by pion and kaon properties and we investigate the quality of generated vector meson 
masses and decay constants. In Sec. II we outline the framework of the DSE approach we 
employ along with the truncation and the Ansatz we use to specify the kernel (or effective 
gluon 2-point function) for both the quark DSE and the bound state BSE. Our investigations 
are conducted with a variation of the kernel Ansatz that was developed in Ref. [jlO| for the 
pion and kaon. To facilitate the analysis and solution of the vector BSE, we have employed 
a convenient set of eight Dirac covariants that satisfy both the CPT constraints and a trace- 
orthogonality property. These are presented and discussed in Sec. III. Also in that Section 
we outline the technique of expansion of the amplitudes in terms of Chebyshev polynomials 
that is sometimes used to resolve the angle dependence and reduce the BSE to a set of 
one-dimensional equations. In Sec. IV the meson decay constants treated here are defined. 
Results are presented and discussed in Section V, and a summary and conclusion follows in 
Sec. VI. Some technical details are collected into an Appendix. 



II. DYSON SCHWINGER EQUATIONS 

In a Euclidean space formulation, with {7^^,71/} = 25^,^, 7^^ = 7^ and a ■ b = J2t=i '^i^ii 
the DSE for the renormalized dressed-quark propagator is 



S{p) 



-1 



Z,i^-p + Z,m{ii) + Z,j g^D^,ip-q)^^^S{q)T:{q,p), (1) 



where D^j^^ik) is the renormalized dressed-gluon propagator, r"(g;p) is the renormal- 
ized dressed-quark-gluon vertex, and = f^d^q/{2TT)'^ represents mnemonically a 
translationally-invariant regularization of the integral, with A the regularization mass-scale. 
The final stage of any calculation is to remove the regularization by taking the limit A — > 00. 
The solution of Eq. (Q) has the general form 

Sip)-'=t^-pAip\fx') + Bip^f^'), (2) 

and the renormalization condition is 

S{py^ = i-fp + m{fi) , (3) 

at a sufficiently large spacelike /i^, with m{n) the renormalized quark mass at the scale /i. 
The renormalization constants for the quark-gluon-vertex, the quark wave-function, and the 
mass, namely A^), Z2(/i^, A^) and Z4(/i^, A^) respectively, depend on the renormaliza- 

tion point and the regularization mass-scale. In Eq. (|I]), S, T"^ and m(/i) depend on the quark 
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flavor, although we have not indicated this explicitly. However, in our analysis we assume, 
and employ, a flavor-independent renormalization scheme and hence all the renormalization 
constants are flavor-independent. 



A. Meson Bethe Salpeter equation 

The renormalized, homogeneous BSE for a bound state of a quark of flavor a and an 
antiquark of flavor h having total momentum P is given by 

^tAv. P) = f^4^(P^ ^' + VP)rfs{q; P)S\q - vP) , (4) 

where rj + fj = 1 describes momentum sharing, T'^j{p; P) is the BS amplitude, and M spec- 
ifles the meson type: pseudoscalar, vector, axial-vector, or scalar. In this paper we consider 
the pseudoscalar and vector amplitudes only. The kernel K operates in the direct product 
space of color and Dirac spin for the quark and antiquark and is the renormalized, ampu- 
tated qq scattering kernel that is irreducible with respect to a pair of qq lines. It is often 
convenient to express Eq. in the abbreviated form 

[TtiP;P)],^ = r Ki:ip,q;P) [xtiQ-.P)],, ' (5) 

Jq 

where Xm^I'^P) •= S"'{q+)^'m{Q'i P)S''{q-) is the BS wave function, g+ := q + rjP, g_ : = 
q — f]P, and the labels r,. . . ,u represent color- and Dirac-matrix indices. This equation 
defines an eigenvalue problem with physical solutions at the mass-shell points P^ = —rri^ 
with m being the bound state mass. 

The canonical normalization condition of the solution of the homogeneous BSE is 

= ^{ TrcD -K) S\q + r/P) ^'(g; K) S\q - r^P)] + 

[x'S,{t,-K)\^^KZ{k,q;P) [xt{T,K 



A r\ 

(6) 



q J k 



where f ^/(fc, — P)* = C~^TM{—k, —P)C, in which C = 7274 is the charge conjugation 
matrix, and X* denotes the matrix transpose of X. The trace in the first term is over both 
color and Dirac indices. If the quark- antiquark scattering kernel K is independent of the 
total momentum P, as is the case in the ladder truncation we consider here, then the second 
term vanishes. 

B. Ladder-rainbow truncation 

We use a ladder truncation for the BSE 

Ki:{p,q;P)^-g{{p-qr)D^-{p-q)lYlA ^lyT.j , (7) 
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which is consistent with a rainbow truncation for the quark DSE 



^1 g'D^UP - q)Y^,Siq)Kiq,p) g{{p - qf) D^-^ip - g)y7M5(?)y7. . (8) 

Here D^^l'^{k) is the perturbative gluon propagator in Landau gauge. The model is completely 
specified once a form is chosen for the "effective coupling" Q{k'^). 

The consistency of Eqs. (|^ and (P) lies in the fact that the axial- vector Ward-Takahashi 
identity is preserved |p],p!0|. This ensures that in the chiral limit the ground state PS mesons 
are massless even though the quark mass functions are strongly enhanced in the infrared. 
In the physical case of explicit chiral symmetry breaking, it also ensures an exact relation 
between the PS meson mass and weak decay constant, the current quark masses, and the 
residue at the PS meson pole in the PS vertex |l|,p!0|. The analysis in Ref. Jl^ shows 
that the next-order contributions to the kernel in a quark-gluon skeleton graph expansion, 
have a significant amount of cancellation between repulsive and attractive corrections for 
pseudoscalar mesons. Indications are that this is also the case in the vector channel, which 
strongly supports the use of ladder truncations in these cases. 

In choosing a form for Q{k'^) we know that the behavior of the QCD running coupling 
a{k'^) in the ultraviolet, i.e. for k'^ > 2-3 GeV^, is well described by perturbation theory. 
In principle, constraints on the infrared form of Qik"^) can be sought from studies of the 
DSEs satisfied by the dressed gluon propagator, D^i,{k), and the dressed gluon-quark ver- 
tex r"(g,p). The latter is often represented by an Ansatz; there is almost no information 
available from DSE studies; the gluon propagator has been often studied via its DSE. If 
the ghost loop and the quark loop in the gluon DSE are unimportant, then the qualita- 
tive conclusion from such studies is that the gluon propagator is significantly enhanced in 
the infrared and well-represented by an integrable singularity such as a regularization of 
1/k'^ Ijl^. Phenomenological studies containing such an enhancement show that dynamical 
chiral symmetry breaking and quark confinement follow in a straightforward and natural 
way from the quark DSE with an empirically correct value for the chiral condensate (qq)^ 
and an excellent description of pion and kaon properties |jlO| . 



Recent gluon DSE studies that include the ghost loop but not the quark loop have 
suggested a weak infrared strength that vanishes at fc^ = for the transverse component of 
Dfj,u{k) due to a strong infrared enhancement of the ghost propagator P^JTPj] . In some studies 



of this type, unphysical particle-like singularities occur in the Ansatz for the dressed ghost- 
gluon and 3-gluon vertices [|T^]. It is apparent that such gluon DSE studies are presently 
limited by the type of truncation that can be accommodated and the preliminary nature 
of the Ansatze employed for some of the dressed vertices. Several lattice studies of D^^{k) 
have been interpreted in terms of an infrared behavior less singular than l/k"^ [20]. The 



phenomenological implications of either type of non-singular infrared behavior for D^^{k) 
have recently been explored within the quark DSE |21[]. It was found that dynamical chiral 



symmetry breaking as represented by a nonzero chiral condensate is either absent or is a 
small fraction of what is required to explain pion phenomena; the produced quark propagator 
does not show quark confinement. 

To provide a quark DSE-based description of pion and kaon phenomena as a basis for 
exploring vector meson properties, we utilize a variation of the following Ansatz introduced 



in Ref. 10 
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fc2 



47r2 



+ An ■ 



r+(l + A:VA2^CD)' 



(9) 



with J^{k'^) = [1 - exp(-fcV[4m2])]/P, r = - 1, and 7™ = 12/(33 - 2Nf). This Ansatz 
preserves the one-loop renormahzation group behavior of QCD for solutions of the quark 
DSE. In particular, the correct one-loop QCD anomalous dimension of the quark mass func- 
tion M(p2) is preserved in its ultraviolet behavior for both the chiral limit {m{fi) = 0, 
anomalous dimension 1 — 7^) and explicit chirally broken case {m{fi) = 0, anomalous di- 
mension 7m). This asymptotic behavior, a characteristic of QCD, is confirmed by analysis 
of the numerical solution in the ultraviolet as described in detail in Ref. |jlO|- The main 
qualitative feature of Eq. is that the phenomenologically required strong infrared en- 
hancement in the region - 0.5 GeV^ is distributed over an integrable S'^^k) singularity ||2^ 



and a finite-width approximation to 6'^{k) normalized so that both terms have the same 
J d'^k. The last term in Eq. (^) is proportional to a{k'^)/k'^ at large spacelike k"^ and has no 
singularity on the real k"^ axis. The parameters u and were not varied freely in the study 
of Ref. |T0|; the fixed values = 0.5 GeV and u = 0.3 GeV were chosen mainly to ensure 
that Qik"^) ~ 47ra(A;2) for k"^ > 2GeV2. The free parameters were D and the renormalized 
u/d- and s-quark current masses to obtain a good description of tt and K properties. 

For the present study of vector mesons, we eliminate the (5-function term from Eq. 
and allow the second (finite-width) term to carry all of the infrared strength. Solutions of 
the rainbow DSE for the quark propagator, when investigated, usually reveal a non-analytic 
behavior in the complex p'^-plane often in the form of complex conjugate branch points p3|j2^ 
that are modified or even eliminated when the gluon-quark vertex is dressed p5|. Subsequent 
use of the propagator solutions in the BSE for the bound state meson should be accompanied 
by a determination that such non-analytic points (that are likely artifacts of the truncation) 
lie outside the complex domain of integration that naturally arises in the search for a solution 
of the BSE in Euclidean metric. The mass of the meson determines the extent of the required 
departures from the quark real axis and the pion and kaon solutions from the Ansatz of 
Eq. (|^) are free of such problems. However, with the parameters of Ref. |]TD[, we have 
found this not to be the case for the more massive vector solutions. The removal of the 
5-function term allows parameters to be easily found to preserve the quality of the pion and 
kaon description while allowing numerically accurate BSE solutions for the vector masses 
reported here. 

We therefore employ the Ansatz 



k^ 



Dk'^e 



+ Air 



iln 



r+{i + kyAic^y 



(10) 



As 



m 

^^QCD 



the earlier pion and kaon studies, we use mt = 0.5 GeV, r = e — 1, Nf = A, 
= 0.234 GeV, and a renormahzation point fi = 19 GeV, which is sufficiently per- 
turbative to allow the one-loop asymptotic behavior of the quark propagator to be used 
as a check. We consider three parameter sets characterized by three different values of u. 
For each parameter set, D is treated as a phenomenological parameter, which was fitted, 
along with the renormalized current quark masses, to obtain a good description of m.,r/K 
and /tt. Subsequently, the vector meson sector was studied without parameter adjustment. 
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For comparison we also report, for the Ansatz of Ref. [|T0|, vector meson masses estimated 
by an extrapolation of the BSE eigenvalue to the mass-shell point. 



III. VECTOR MESON BETHE SALPETER AMPLITUDES 

The general form of a vector vertex r^(g; P) can be expressed as a decomposition into 
twelve independent Lorentz covariants, made from the three vectors 7^, the relative momen- 
tum q^, and the meson total momentum P^, each multiplied by one of the four independent 
matrices 1, 7 ■ g, 7 ■ P, and a^^q^P^. Since a vector meson BS amplitude is transverse the 
number of allowed covariants reduces to eight, so that the general decomposition of the 
vector BS amplitude is 

8 

rj:(9; P) = Y1 P) ^^(^'' ^ ■ ^' P") ' (11) 

i=l 

with the invariant amplitudes Fi{q^, q ■ P; P^) being Lorentz scalar functions. The choice for 
the covariants T^{q; P) to be used as a basis is constrained by the required properties under 
Lorentz and parity transformations, but is not unique. The BSE Eq. (|) must be projected 
onto the covariant basis to produce a coupled set of eight linear equations for the invariant 
amplitudes Fi to be cast in matrix form. This requires a procedure to project out a single 
amplitude from the general form Eq. (0). It is therefore helpful if the chosen covariants 



satisfy a Dirac-trace orthonormality property. 

We have chosen the following set of dimensionless orthogonal covariants 

Tl{q-,P)=ll, (12) 

TlW, ^) = ^ ■ <i) - bli'i'f) > (13) 

n{q:P) = ^[ql{l-P)) , (14) 

n{q; P) = '-^ [llil ■ P){1^ ■ q) + ^ (7 ■ P)) , (15) 
T'^{q;P) = -^ql, (16) 

^ (?; P) = ^ Kii^ ■ <i) - (7^ ■ > (17) 

q V Z 

Tliq; P) = (1 - cos^ e) (7j (7 ■ P) - (7 ■ PH) - ^ T'M, P) , (18) 

T,%P) = ^^gJ(7^-g)(7-P), (19) 

where V'^ is the component of V transverse to P 

^; = ^.-^^^^, (20) 
8 



and q ■ P = q P cos 9. Note that at the mass-shell P = im. The orthonormality property 
satisfied by these covariants is 

^Trc P) r^(g; P)] = ^cos 9)5., , (21) 

where the functions fi{z) are given by fi{z) = 1, fi{z) = |(1 — z^) for i = 3,4,5,6 and 
fi{z) = |(1 — 2;^)^ for z = 2, 7, 8. For later use we also note the relation 

TT 

d9 8111^9 fi{cos 9) = 1- (22) 

The covariants are dimensionless and independent of the magnitudes q and P. These prop- 
erties are helpful in allowing the relative magnitude of the amplitudes to be a qualitative 
measure of the dynamical importance of the various covariants. A more quantitative mea- 
sure can depend on the particular observable being studied; amplitudes that are unimportant 
at low momenta can become dominant when high momentum behavior of the bound state 
solution is being probed. 

For unfiavored mesons that are eigenstates of C (charge conjugation), such as the p, uj 
and 0, there is an additional constraint on the BS amplitude^ to obtain a specified C-parity. 
Of the eight covariants given in Eqs. (|12[)-(|19]), and are even under C, the others are 
odd under C. The only remaining quantity that can produce a desired uniform C-parity is 
q-P which is odd under C . Thus a C = — solution (such as the p and (p) will have amplitudes 
F3 and that are odd in g ■ P while the remaining amplitudes are even in g ■ P. For the 
fiavored vector meson K* , which is not an eigenstate of C, each amplitude will contain both 
even and odd terms in g-P. Since the ladder truncation of the BSE is invariant under charge 
conjugation if equal momentum sharing (77 = 0.5) is used, the observation of the above odd- 
even behavior in g • P of Pj can be used as a test of numerical accuracy. Alternatively, the 
amplitudes Pj can be expanded in terms of a basis of functions that are appropriately odd 
or even in cos 9 to save significantly on computer time and memory. Because the mass-shell 
condition makes the magnitude P imaginary, it is not difficult to verify that with definite 
C-parity, each amplitude associated with our chosen basis of covariants is either purely real 
or purely imaginary. The amplitudes Pj for the K* solution are in general complex due to 
the dependence upon all powers of g ■ P. 

After using the representation Eq. ( pH]) for the solution in terms of the covariant basis, 
followed by projection using the Dirac-trace orthonormality property Eq. (|2l|) , the homo- 
geneous BSE Eq. (^) for a meson with fiavored constituents ah reduces to a set of coupled 
integral equations for the eight functions F°'^{q^, q ■ P; P^) in the form 

Ff{p\p. P; P')f\{z) = 1 Qiip - g)2) D^^^ip - q)Ff{q^ q ■ P; P') x 

^Trz5 [r;(p; Ph.S^iq + r/P)T^(g; P)S\q - r/P)7.] . (23) 



^We do not discriminate between up and down quarks, and do not take into account electromag- 
netic corrections; therefore the BS amplitudes for are equal to those for p^. Furthermore, the 
ladder truncation cannot discriminate between isovector and isoscalar mesons; therefore the p and 
the Lo are degenerate in this truncation. 



9 



The above system of equations was solved by two complementary methods. The first method 
was a direct treatment as an integral eigenvalue equation A(P^)F = K(P^)F for a set of 
functions F of two variables: and z = cos 6. An iterative method is used to determine the 
smallest m satisfying A(— m^) = 1. Both variables were discretized via Gaussian quadrature 
and the summations for the double integration were carried out at each iteration. This has 
a high demand on computer memory. 

In the second method, the angle dependence of the amplitudes is expanded in the form 

oo 

q.p-P^) = Y^ W,{q'; P') f/,(cos0) , (24) 

where the Uj{z) are Chebyshev polynomials of the second kind. This allows the angle 
integrations in Eq. ( p3l) to be carried out to produce an integral equation in one variable but 
for a larger set of functions ^Fi{q^] P^). For C = — eigenstates such as p and amplitudes 
P3 and Fq will require only odd order Chebyshev terms while the other amplitudes will 
require only even terms. In practice, the number of Chebyshev terms required is quite low 
(one or two terms) so that the memory requirements are effectively reduced in this second 
method. The solutions from the direct two-variable approach can be projected onto the 
Chebyshev basis as a check on the second method and also as a means of presentation. 

The specific normalization condition for the vector meson solutions of the ladder BSE 
follows from Eq. (^ and is 



d N 

= ^ Y Trn R"(g; -K) S'^{q + ^P) (g; K) S\q - r/P)] 
where the factor 1/3 appears because the three transverse directions are summed 



, (25) 



IV. ELECTROWEAK DECAY 



Here we summarize the definition of, and our convention for, the vector meson leptonic 
and electromagnetic decay constants and their explicit relationship to the BS amplitudes. 
The electromagnetic decay mediated by a photon (e.g. p°, cj, 0), and the leptonic decay 
mediated by a W-boson (e.g. p^, K*^)^ are described by the vector decay constant defined 
by 12§ 

hmyef{P) = {0\(fj,q'^\V'^\P,X)) , (26) 

where e|t^^ is the polarization vector of the vector meson satisfying e*^^^ ■ P = and normalized 
such that e^^''* e^i"^ = 3. This is completely analogous to the definition 

fpP, = (0|g%^75g1P"''(P)) (27) 

for the pseudoscalar decay constant that corresponds to fj^ = 131 MeV under the normal- 
ization convention of Eq. (^). The vector decay constant from Eq. ( PBj ) can be expressed as 
the loop-integral 
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7 N H'^n 

fvmy = ^J Trz, [i,S'^{q + vP^iq; P)S\q - vP)] , (28) 

which is exact if the dressed quark propagators and the meson BS amphtude are exact 0. 
In the next section we use Eq. (^) to calculate the decay constants fp, and fx*- 
The coupling of the 

|p°) = i= {\uu) - \dd)) (29) 

to the photon is conventionally expressed via a dimensionless coupling constant gp in the 
form 

2 

771 _ 

^6W(P) = (0|QQ7,Q|P°(P,A)), (30) 
9p 

where the flavor multiplet of quark field spinors is Q = column(M, d), and Q is the quark 
electromagnetic charge operator. The normalization condition given in Eq. ( p5l) is in a form 
appropriate for a single flavor configuration QaQb, not for a multi-flavor configuration state 
like the p°. For such states, Eq. (^) can be generalized by promoting the quark propagators 
to flavor matrices S = diag(5'", S'^), multiplying BS amplitudes Tp by the appropriate flavor 
matrix, and tracing over flavor indices as well. The isospin-symmetric limit with 5" = S'^ 
produces BS amplitudes that are independent of flavor labels; the p'^ BS amplitude, for 
example, can then be expressed as (rs/v^) E^ where E^ is the normalized BS amplitude for 
the p^. Use of Eq. (|2|) in Eq. (H) then gives 

Tr^ [^pS^=\q + vP)^f{q; P)S^=\q - r^P)] 

(31) 



ml Z2N, r"- dS 



gp 3V2 J (27r) 
_ fp^p 



V2 • 

The decay width TpO^e+e- = 6.77 keV leads via 



An Op' 77lr. 

~^p 



'±/i LX nip /on\ 
FpO^e+e- = (32) 



to the value Qp = 5.03, that is fp = 216 MeV. Note that the isoscalar version of these con- 
siderations produces an extra factor of 1/3 on the right of Eq. (|3lD for the coupling of the 
a; to a photon. The partial width F^^^e+e- is indeed about 10 times smaller than FpO^e+e-- 
In a similar way, the coupling of the photon to the 0, assumed to be a pure ss-state, is 
defined as 



2 
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^eW(^) = |(ON-7.^|0(P,A)), (33) 
9(j) 



and the relation between and the vector decay constant is 
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9<p 



Tin [i,S'{q + vPK'S\q-r]P)] 



(34) 
(35) 



The partial width of the 



e~^e decay is 



(36) 



and the experimental value 1.37 ± 0.05 keV produces = 237 MeV, that is = 12.9. 

The decay constant fv determines not only the coupling of the neutral vector mesons to 
a photon, but also the coupling of and K*^ to the weak vector bosons W^. There are 
no data available for the leptonic decay of these charged vector mesons, but the couplings 
can be extracted indirectly from the decays t ^ pvr and r — > K* Vr. The partial width for 
such a decay is 



mt 



1 + 



2 ml 



(37) 



With the experimental values for the partial decay width |^ ^r^pv^ = 25.02% Ftotai and 
Tr^K*vr = 1-28% Ttotai, and the CKM matrix elements = 0.974 and Ks = 0.220, this 
gives a ratio 



K* 



1.042 



(38) 



and thus a decay constant fx* = 225 MeV, if we use the experimental value /„ = 216 MeV. 



With the available data, the absolute value of fp using Eq. ( pTD gives fp = 208 MeV. We 
expect however that the direct determination of fp through —>■ e+e~, giving fp = 216 MeV, 
is a more accurate determination of this decay constant. In particular, most higher-order 
corrections to the electroweak vertex are likely to cancel in the ratio of the partial decay 
widths, and therefore we use the ratio in Eq. (p8|) to extract the experimental fx*- 



V. NUMERICAL RESULTS 

In Fig. ID we show our Ansatz for the effective interaction, Eq. ([10|) , for the three different 
parameter sets we have explored, characterized by the values of ut, together with the 1-loop 
perturbative coupling for comparison. We use three different values of the parameter u, 
constrained only by the requirement that the perturbative coupling above = 3 GeV^ 
should be reproduced. It is only in the infrared region, below q^ = 2 GeV^, that there is a 
significant difference between the three parameterizations and the perturbative result. The 
parameter D and the current quark mass m^/dip) are fixed by fitting m^^ and /jr. Next, the 
strange quark mass ms{p) is determined by a fit to the kaon mass. The resulting value of the 
kaon decay constant fx is within 3% percent of the experimental value, almost independent 
of the parameter set for the effective interaction. All three parameter sets lead to a good 
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description of the pion and kaon masses and decay constants, as well as a reasonable value 
of the chiral condensate. In Table | we have summarized these results for the three different 
parameter sets, together with the results from Ref. [p!0|| . 

With our parameterization, the quark mass function M(p^) = B{p'^) / Aijp') has qualita- 



tively the same behavior as obtained in Ref. |10|. With a Euclidean constituent-quark mass 



M defined as the solution of p = M (p ), we obtain constituent quark masses of about 
M^/d = 300 - 500 MeV for the light quarks, and Ms = 500 - 640 MeV for a strange quark. 



spanned by the three parameter sets; the parameterization of Ref. ||T0[ gives constituent 
masses M^/a = 560 MeV and M, = 700 MeV. 

A. Results for vector meson observables 

In Table |l| we present our results for the vector meson masses and decay constants. 
The full angular dependence was retained in the calculation of these results: we solve the 
set of integral equations Eq. (p3D with the FiijP'^p ■ P; P^) treated as functions of two 
variables p"^ and z = cos 6. This eigenvalue problem defines physical solutions at the mass- 
shell P2 = 

— TTiy. All calculations with the gluon Ansatz of Eq. ( |10| ) were performed at 
the physical mass-shell; the calculations we have performed with the parameterization of 
Ref. [|I^ for comparison involved some extrapolation^ to the mass-shell, which makes these 



results less accurate. In particular the integral for the normalization condition, Eq. (pSf ), is 
very sensitive to such an extrapolation, which is why we do not report the decay constant 
for this particular model. 

All parameterizations we used give equally good results for the masses and decay con- 
stants: the results are fairly insensitive to changes in uj and D, as long as they are fit to ttIt^, 
fir and rriK- Our result for nip is typically 5% too low, whereas mx* and are typically 
5% too large. Our result for the decay constants are within 10% of the experimental value 
for fp and fx*, and within 10% to 15% for f^, depending on the parameter set. This agree- 
ment with experiment is quite encouraging, given the fact that the parameters are fixed by 
pseudoscalar observables. 

From Table ^ we can also conclude that only five of the eight covariants are qualitatively 
and quantitatively important for the vector meson masses and decay constants; this seems 
to be general, i.e. independent of the parameter set used. Of course, the relative importance 
of different covariants in a BS amplitude does depend on the observable under consideration. 
Also, use of a basis set of eight independent covariants that is different from the present 
basis given in Eqs. (p!2D-(p!9D, could produce a different conclusion concerning the number of 
important covariants. 

In Fig. we show the behavior of the leading Chebyshev projection of the invariant 
amplitudes of the p BS amplitude, ^F^{q^; P^). This and the other plots of the BS amplitudes 
are produced with the parameter set u = 0.4 GeV and D = 0.93 GeV^; the results for the 
other parameter sets look qualitatively the same. The leading amplitude for the pion, E^^, 



^The extrapolations were necessary because of nonanalytic behavior of the resulting quark prop- 



agator as discussed in Sec. 11 B 
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and for the rho, Fi, are very similar; however, this similarity might be accidental. Of the 
sub-dominant amplitudes, F4 and F5 are significantly larger than the rest. The magnitude 
of the amplitudes Fq, F7, and Fg is much smaller than that of Fi, F4^, and F5, as is evident 
from Fig. this makes it understandable why these amplitudes contribute so little to the 
vector meson masses and decay constants. From this figure one might conclude that the 
amplitudes F2 and F3 have a similar minor role. However, it turns out that these amplitudes 
are essential for the convergence of the loop integral for the decay constant, Eq. (EBf), as 



discussed below, in Sec. VB 



To study the relevance of the various covariants for physical observables in more detail, 
we calculate the vector meson masses and decay constants using different subsets of the eight 
covariants in our basis. These results are given in Table |T| for one particular parameter set, 
together with the results from use of only the leading Chebyshev moments of each amplitude 
q ■ P; P^). Note that the leading Chebyshev order for the K* is zeroth order for all 
amplitudes Fi, in contrast to the case for the p and 0: the functions F3 and Fq are odd in 
q-P for the p and because of charge conjugation symmetry, so the leading Chebyshev order 
is Ui{cos6) for those mesons. It is evident that for the p and (p only the leading Chebyshev 
moment is needed to get accurate results for the masses and decay constants; but the second 
Chebyshev moment of Fi is needed for strict convergence of Eq. (PB]). We expect this to be a 
general phenomenon: practical calculations of hadron observables might be facilitated by a 
suitable parameterization of the leading Chebyshev moments of the amplitudes Fi through 
F5. For the K*, which is not a charge conjugation eigenstate, one needs at least the zeroth 
and the first Chebyshev moments for an accurate description. 

Another difference between the K* and the p and mesons, is the dependence on the 
momentum sharing parameter 77 in Eq. (|2^) . Charge conjugation dictates use of = 0.5 
for the p and the 0. For the K* there is no such constraint and we explored momentum 
partition sets {rjuifis) varying between (0.5,0.5) and (0.4,0.6). Physical observables are in 
principle independent of this partitioning; any dependence of K* physical observables on 
{Vu,Vs) would signal an inadequacy of the ladder truncation or subsequent approximations. 
We find that the results for mx* and fx* are indeed unchanged under variation of the 
momentum sharing, as long as all covariants and the full angular dependence are taken 
into account. Once certain amplitudes are dropped and/or the angular dependence of the 
amplitudes is truncated, physical observables do become dependent on {riu,fjs)'- variations 
between (0.5, 0.5) and (0.4, 0.6) lead to changes in mx* and fx* of up to 5%. 

A comparison of the BS amplitudes of the three different vector mesons is made in Fig. |^. 
This figure clearly shows the difference between the p and mesons on the one hand, and 
the on the other: while the leading Chebyshev moments of the p and amplitudes are 
very similar to each other and to the corresponding moments of the K* amplitude, the latter 
has both even and odd moments, due to the lack of C-parity. This is especially evident for 
the amplitudes F^{q, q ■ P; P^) and F^lq, q ■ P; P^), which have no zeroth Chebyshev moment 
in the case of the p and 0, but have a significant zeroth Chebyshev moment for the K*. 

B. Asymptotic behavior of the BS amplitudes 

The asymptotic behavior of the BS amplitudes for the p-meson is shown in Fig. ^, and as 
in the pseudoscalar case, all amplitudes behave like or 1/q^, up to calculable logarithmic 
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corrections. We emphasize that in QCD these logarithmic corrections are essential for the 
convergence of the integral for the decay constant. Evaluation of the trace in Eq. (^) for 
equivalent flavors and equal momentum partitioning gives the leading behavior 



fvmv 



ZoN, 



3 

- 32 q' (1 - cos^ efa+a- F^i- . .)/V5 - 16 q' cosO (l - cos^ 9) a+a^, F^{. . .) 
+ t8V2 Pq (1 - cos^ e) (T>; F4(. . .) - ^8 g (l - cos^ 6) + F,{. . .) 

+ O ((Fe + Fr + Fs)qa^af) } , (39) 

where the vector and scalar components of the quark propagator 

A{q) 

A\q)q^ + B\q)' 
B{q) 

~ A\q) q' + B^q) ' 

and := f\q±)- The last terms in Eq. (|5PD, proportional to Fe, Fy, and Fg, give small and 
convergent contributions, since they behave in the ultraviolet as 



FM\q-P-,P^)qoU^ ~ -, (42) 

up to logarithmic corrections. Both contributions involving F4 and F5, which fall off like 
1/g^, are also ultraviolet finite. Since the amplitudes Fi, F2, and F3 fall off as 1/g^, simple 
power counting shows that their individual contributions to Eq. (^) are logarithmically 
divergent, even accounting for the cut-off dependence of Z2(A^,yU^)2. 

In order to analyze the asymptotic behavior produced by the vector meson BSE in 
more detail, we follow the strategy used in Ref. for the asymptotic behavior of the 
function -B(p^) from the quark DSE. The key step is to replace the effective running coupling 
Q{{p — qY) by ^(max (p^, g^)). In the ultraviolet region the running coupling behaves like 
l/ln(?/) with y = /e^/Aq^^j^ which is a slowly varying function; therefore the error made in 
using this approximation is under control. In the infrared region, such an approximation is 
not to be trusted. After this approximation, and with use of the Chebyshev decomposition 
for the angular dependence of Fj(g^,g ■ P,P') in Eq. (^31) , all angular integrations can be 
performed analytically. For the leading Chebyshev moments of the p BS amplitudes, Eq. (^3]) 
produces integral equations of the form 

F^{x) = ^ r dy¥:^^^{x,y)F,{y) + ^^ H dy Kl^Ux,y)^ , (43) 
m X Jq m y 

where x = pVAqcd> V = iV^qcb^ H^) = %{q';P') for i = 1,2,4,5,7,8 and Fi{x) = 
^Fi{q']P'') for i = 3,6. Now the coupled integral equations can be converted to a set 



^The factor Z2{A'^,ij?) ensures gauge invariance and cancels logarithmic divergences in covariant 
gauges other than Landau gauge. 
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of coupled linear differential equations, which can be solved in the ultraviolet region by 
assuming a series expansion in both x and Inx. For completeness, we have given the relevant 
kernels K*-' and other details in Appendix 

The analysis in Appendix ^ shows that the ultraviolet behavior of the amplitudes Fi (x) , 
F2{x) and F^^x), is of the form 

^fl + f:c,(ln.)-A, (44) 



Fix, 

X 



and the steps leading to identification of the power a and the leading coefficients Oj are also 
given there. The leading ultraviolet behavior is found to be 

'F,{q'-P') = (45) 

= i^.(x)-^^^^|f^, (46) 

%{q';P') = F3(a;)~^i^, (47) 

with 

a = -l + 7^/108 . (48) 

The overall constant ai is not determined by the homogeneous BSE; its value follows from 
the normalization condition. 

Our numerical results show that the leading ultraviolet behavior of the BS ampli- 
tudes is governed not only by the leading Chebyshev moments ^Fiiq"^; P^), °F2(?^; P"^), and 
^F^iq'^] P^), but also by the second Chebyshev moment ^Pi(g^; P^), see Fig. ^ Numerically, 
we find in the ultraviolet 



'F2{q';P') 

0Pl(g2;P2) 

'Fs{q';P') 

0Fl(g2;P2) 

'F,{q';P') 



0.48 ±0.01, (49) 
0.33 ±0.01, (50) 
-0.11 ±0.005, (51) 



while all other ^Pj fall off faster. This is in excellent agreement with the analytical results 
for the relative magnitudes of the leading Chebyshev components, Eqs. (|45| ) -(|47| ) . The 
power a of the logarithm is much harder to determine numerically; our results indicate 
—0.95 < a < —1.0, which is consistent with a = —0.996 from Eq. (^). We have not studied 
whether the inclusion of ^Pi in the analysis of the asymptotic behavior would change our 
analytical result for a; our numerical results indicate that it will not influence the coefficients 
tti of °Pi, °P2, and ^Ps, nor will it change the power /3 in Eq. (|A7|) . 

The ultraviolet behavior of the integral for the decay constant, Eq. (P5|), can now be 
analyzed in more detail. The ultraviolet behavior of the functions Pi, P2, and P3 does indeed 
lead to individual divergent integrals for a > —1. However, the combined contribution is 
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fvmvr^J rfy^^^ |6ai-4v^a2-4a3 + 26i| , (52) 

where bi is the coefficient of the second Chebyshev moment ^Fi(g^;P^), that is, the coun- 
terpart of ai in Eq. (^51) . Use of the asymptotic behavior we have found analytically, 
02 = ai2\/5/9 and 03 = ai/3, shows that the integral for the decay constant is finite if 
bi = — ai/9, which agrees with our numerical result, Eq. (|5T|). This cancellation between 
naive divergences coming from different covariants provides an illustration of how renor- 
malizibility is realized; it is expected since the one-loop renormalization group behavior of 
QCD is preserved in our rainbow-ladder truncation of the DSE and BSE. It is the vector 
counterpart of a similar cancellation in the integral for the pseudoscalar decay constant: 
numerically [|1^, it was found that in the ultraviolet region the tt BS amplitudes satisfy 



= 2F^/g^, which makes the integral for finite, although the separate contributions 
from and diverge. The above analysis, when applied to the pseudoscalar BSE, pro- 
duces an asymptotic behavior of the amplitudes F^^ and that exactly gives = 2-F^/p^. 

Finally, in Appendix 0, the same ultraviolet analysis applied to the vector BS amplitudes 
F4 and F5 gives 

F.{x)- . , (53) 



^3/2 ^Iyix 



In principle, the influence of Fq, F7, and Fg might change the power of the logarithm for F4 
or F5, but we expect no change in the leading ultraviolet behavior of Fi, F2, and F3. 



VI. CONCLUDING REMARKS 

We have calculated the light vector meson masses and the decay constants associated 
with electromagnetic and leptonic decays using the ladder truncation for the meson BSE 
in conjunction with the rainbow truncation for the quark DSE. We use an effective quark- 
antiquark interaction Q{k'^)/k'^ with one phenomenological parameter, which is fitted to 
reproduce the two other parameters are the current quark masses niu/d and rus which 
are fixed through and m^- The calculated values for the vector meson masses are within 
5% of the experimental values; the decay constants are within 10% of their experimental 
values. These results are fairly robust: they are weakly dependent upon the scale at which 
the interaction starts to deviate from the perturbative behavior, as long as the parameters 
are fitted to pseudoscalar observables. 

An earlier BSE study in a related framework produced qualitatively comparable 
results for m^, m^* and m,^ as part of a study that included heavy mesons and incorporated 
five quark flavors. Vector meson decay constants were not considered. That approach 
produced a dependence upon the momentum-sharing parameter 77 that is stronger than 
what we find. The present results for physical observables, such as the mass and decay 
constant, are independent of the momentum sharing, as long as all relevant covariants and 



the full angular dependence are included in the calculation. A recent work |29| has explored 



the feasibility of extracting ground state vector meson masses from the large Euclidean time 
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behavior of the quark current- current correlator as calculated from the ladder truncation of 
the inhomogeneous BSE for the vector vertex. Only the p was studied, the BS amplitudes 
were not extracted and the decay constant was not calculated. 

Of the eight allowed transverse covariants, five are both quantitatively and qualitatively 
important, whereas the remaining three amplitudes contribute little to the mass and decay 
constant. Neglect of these three amphtudes changes the calculated masses by only 2% and 
decay constants by 8%. For the p and (f) the dependence of the BS amplitudes on g ■ P 
is very small; truncation to the leading Chebyshev moments leads to very similar results. 
However, the second moment of Fi is needed for convergence of the loop integral for the 
decay constant. This suggests that, in general, hadronic observables can be well-described 
by a rather limited number of covariants and Chebyshev moments. For the however, 
more Chebyshev moments are required, since it is not a charge conjugation eigenstate. Our 
numerical results can be used to guide the development of approximating forms for the BS 
amplitudes for calculation of a variety of observables such as electromagnetic form factors 
and strong decays typified by as p — tttt and KK. 

The ladder truncation of the BSE is known to be a good approximation for flavor nonsin- 
glet pseudoscalar mesons [|17|, and it is expected to be reliable for vector mesons as well. This 
is to be contrasted with the scalar channel, where the same analysis revealed |^D|] that the 
next-order corrections are much more important. For fiavor-singlet mesons, there are also 
contributions from diagrams corresponding to quark annihilation to time-like gluons. These 
play an important role for pseudoscalars, e.g. in the generation of the rj' mass through the 
axial anomaly [0]. For the fiavor-singlet vector mesons however, there is no such anomaly. 
Also, if quark annihilation diagrams play a major role for such mesons one would expect 
more flavor mixing than is evident for the uj and 0. It is therefore reasonable to expect the 
ladder truncation to be appropriate for vector mesons. For the ground state vector mesons 
considered here, there is an open decay channel to a pair of pseudoscalars (e.g. p tttt), 
but this is a P-wave coupling that tends to suppress the mechanism relative to such a decay 
of a scalar. Estimates of the effects of meson loops on the p mass vary between 2% and 



10% With the BS amplitudes calculated here we expect to be able to investigate the 
effects of meson dressing more accurately in the future. Note that both the meson dressing 
and the quark annihilation diagrams can contribute to the splitting between the p and uj, 
which are degenerate in the ladder truncation. 

The task of modeling vector mesons within QCD at finite temperature and chemical 



potential has recently begun with extremely simplified Ansatze for the kernel of the BSE 



The present work may provide valuable guidance for the extension and improvement of 
such efforts to explore the behavior of vector qq states and correlations relevant to chiral 
restoration and quark deconfinement transitions. 
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APPENDIX A: ASYMPTOTIC BEHAVIOR OF THE BS AMPLITUDES 



In order to analyze the asymptotic behavior of the BS amphtudes '^Fj{q^\ P^), we have 
to perform all angular integrals analytically. These angular integrals have measure 



dOp sin 6p 



dOq sin 6q 



2tt 



sm I 



(Al) 



where 6p is the angle between the external momentum p and P, and 6g is the angle between 
the integration momentum q and P. This type of integral can be performed with the help 
of the appendix of Ref. ||33|, and some typical results are 

p ■ q pq min(p, q) 



dn 
dn 



{p — qy 2max(p, g)3 ' 
p ■ q pq min(p, q) 



ip ~ max(p, g)3 (max(p, qY — min(p, qY) 



(A2) 
(A3) 



Other, more complicated, integrals can be expressed in a similar way. A common feature of 
these angular integrals is that they can all be expressed in terms of max(p, q) and min(p, q). 
This allows us to convert the integral equations to differential equations [P5| . 

We have performed all the angular integrals in the five coupled integral equations for 
°Fi(g2;p2^, 0F2(g2;p2), ^F^i^q^ ] P"^) , ^Fi{q^;P'^), and ^F^{q^ ; P'^) , ignoring the functions Fg, 
F7, and Fs, and truncating the Chebyshev moments at the leading order. With the leading 
ultraviolet behavior of the functions Fi, F2, and F3 considered first, the relevant kernels are 



■^a;>j/('^5 y 

^f>y{x^y 
^f>y{x.y 



K22 



x>y 

23 

x>y 



{x,y 



K31 



x>y 



{x,y 



^x>y{-^i y 



1 

4a; ' 



-V5 



6 X 

-1 
6x ' 



-V5{4y-3x) 

5{9y-8x) 
216x2 ' 

-VE(8x-3y) 
216x2 

~{4y-3x) 
36a? ' 

-V5{8x-3y) 
144 x2 

ny-8x 
144x2 5 



mxi^.y 
^y>x{-^i y 

^y>x{-^i y 
^y>x{-^i y 
^y>x{-^i y 
^y>x{-^i y 
^y>x{-^i y 



These are to be inserted into the integral equation Eq. 



X 

-V^jiy-x) 
18 y2 

-(4 y-x) 
18?/2 5 

-VE'. 



54 y^ ' 

5x2 
216 ' 



216 y^ ' 

— X 



36 s/2 ) 



IMy'-^ ' 

(25x-16y)x 
144 y3 

3f), which is 



Fix) 



In 



Inx 



^x>y 

(x, y)Fj{y) +'ym dy K^>j.(x, y) 



In?/ 



(A4) 



(A5) 



where Fi{x) = ^Fi{q'^;P^) for i = 1,2 and F3(x) = ^F^{q^] P'^). This set of coupled integral 
equations can now be converted into a set of coupled fourth-order differential equations for 
Fi_3(x) of the type 
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X^K4i-Fj {x)+X^KsiF^ (x) + x'^ K2iF^ (x) + X KiiFj^{x) + KoiF{x)i = . (A6) 

Substitution of the series expansion 

F.{x) = + ^^"'^^") ^^^^ 

into the set of differential equations leads to a set of coupled equations for the powers a and 
/3, and the leading coefficients Oj. It is easy to see that all terms in the differential equation 
have the same power of x, and collection of all the leading powers of In x gives an equation for 
the power p. One of the solutions of this equation is P = 1, which is obviously the physical 
solution, see Fig. ^. The next-to-leading order terms lead to three coupled equations for 
the four constants a, ai, 02, and 03; the homogeneous BSE allows for an arbitrary overall 
scaling and we set ai = 1. The solution for the other constants is then 

« = -1 + ^ , (A8) 
108 ' ^ ^ 

a2 = ^ , (A9) 

as = ^ . (AlO) 

Note that the powers a and (3 are the same for all three functions. Differences between the 
functions only arise from differences in the coefficients for the leading, and the subleading 
coefficients c^. 

Next we consider °F4, which decouples from the other amplitudes after performing the 
angular integrals. The only nonzero kernels in Eq. ([A5|) are 



Kiyix,y) = ^, ^Ux,y) = ^. (All) 

The resulting asymptotic behavior can be expressed by Eq. (|A7|) with 

P4 = I (A12) 
a4 = -l + |7m. (A13) 

This is in agreement with the numerical result, see Fig. |. The equation for F5 is more 
complicated, since ^F^ does couple to °Fi, °F2, and ^^3. The relevant kernels are 

K51 „,\ _ M{y) T^5l / \ _ y^Mjy) 

'^^xyyK-^i y ) ~ 2x^/'^ ' '^^y>x\-^iy) ~ 2y^ ' 

K52 „A _ V5M{y) t^52 /„ „,\ _ Vb^Mjy) 

J^x>yl-^) y) - 3^3/2 , '^y>x\-^-: y) - 3y2 , (AIA) 

^53 „,\ _ M{y) T^53 / x _ ^Mjy) 

'^^xyyK-^i y ) ~ 3T.3/2 ) '-\>x\-^^ y J 3„2 ) 



x>yK-^^yj 3a;3/2 ' "-^yyxK-^^yj 3y2 

55 / \ _ yy T^55 ( \ _ 

x>y\-^iy) ~ 2a;3/2 ' ^^y>x\-^-i y ) 2y3/2 • 



However, a careful analysis shows that the leading ultraviolet behavior of ^F^ is not influenced 
by coupling to other amplitudes; the leading behavior arises from K^^ only. The result is 
Eq. (^) with 
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/35 = i, (A15) 
as = -1 + |7m , (A16) 

also in agreement with our numerical results, see Fig. ^. The influence of ^F^ on our previous 
results for °Fi, °F2, and ^F^ can be examined for consistency. Those three amplitudes fall off 
like while contributions from ^F^ to these amplitudes via the differential equations in 
Eq. ( |A6|) will be suppressed by a factor of M{x)/y/x, and thus contribute to the subleading 
behavior only. 
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q' [GeV'] 



FIG. 1. The Ansatz for the effective qq interaction Q{q^)/q^, Eq. (|To|), for the three parameter 
sets, together with the 1-loop perturbative result for comparison. 
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0.5 1 1.5 2 

q [GeV] 

FIG. 2. The leading Chebyshev projections of all eight p BS amplitudes, normalized to 
°Fi(0) = 1, with an effective qq interaction, Eq. (|^), with uj = 0.4 GeV, D = 0.93 GeV^ The most 
important amplitudes, Fi-F^, are labeled by lines with symbols. 
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FIG. 3. Leading and subleading BS amplitudes for the />, iC*, and (f> mesons, (a) zeroth 
Chebyshev projections of -Pi, and for the K* also the first Chebyshev projection, (b) as (a), but 
then for F2, (c) first Chebyshev projections of F3, and for the K* also the zeroth projection, (d) as 
(a), but then for F4, (e) as (a), but then for F5, (f) as (c), but then for Fq, (g) as (a), but then for 
Ft, (h) as (a), but then for Fg. The parameters are the same as in the previous plot: u> = 0.4 GeV, 
D = 0.93GeV2. 
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q' [GeV'] 



FIG. 4. The ultraviolet behavior of the p BS amplitudes: the leading Chebyshev moments of 
Fi-Fs, obtained using the full angular dependence, are shown by the symbols. The lines display 
the analytically calculated behavior for -F1-F5, given by Eqs. (45)-(47), (|5^), and 
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FIG. 5. The ultraviolet behavior of the p BS amplitudes: the ratio ^Fj(g^; P^)/°Fi(g^; P^) for 
the leading amplitudes. 
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TABLES 





experiment 




id — U.O \jG V 


UJ — U.4 (jG V 


UO — U.O \jG V 




(estimates) 


Tiof 111 rfl 
xiei. ||iij| 


U — i.zD \jQ V 


u — u.yo (jG V 


u — u. ( y Lie V 




(0.236GeV)3 


(0.241GeV)3 


0.242 


0.241 


0.243 


">=lGeV 


5 - 10 MeV 


5.5 MeV 


5.54 


5.54 


5.35 


"^Ai=lGeV 


100 - 300 MeV 


130 MeV 


124 


125 


123 




0.1385 GeV 


0.1385 


0.139 


0.138 


0.138 


f. 


0.1307 GeV 


0.1307 


0.130 


0.131 


0.131 


rriK 


0.496 GeV 


0.497 


0.496 


0.497 


0.497 


fK 


0.160 GeV 


0.154 


0.154 


0.155 


0.157 



TABLE I. Calculated values of the properties of light, pseudoscalar mesons, for the param- 
eterization of the effective interaction Eq. (|l^), using three different parameter sets, and also for 
the parameterization of Ref . |l^ . 







P 




K* 












nip 


fp 


niK* 


fK* 




f4> 


experiment 


0.770 


0.216 


0.892 


0.225 


1.020 


0.236 


All amplitudes Fi- 


Fs 




LO = 0.3 GeV, D = 


1.20 GeV^ 


0.747 


0.197 


0.956 


0.246 


1.088 


0.255 


UJ = 0.4 GeV, D = 


0.93 GeV^ 


0.742 


0.207 


0.936 


0.241 


1.072 


0.259 


UJ = 0.5 GeV, D = 


0.79 GeV^ 


0.74 


0.215 


0.94 


0.25 


1.08 


0.266 



amplitudes Fi . . . Fr, only 



Maris-Roberts Ref. 


0.71 




0.95 




1.1 




io = 0.3 GeV, D = 1.20 GeV^ 


0.737 


0.192 


0.942 


0.235 


1.080 


0.247 


UJ = 0.4 GeV, D = 0.93 GeV^ 


0.729 


0.199 


0.919 


0.229 


1.062 


0.250 


UJ = 0.5 GeV, D = 0.79 GeV^ 


0.731 


0.207 


0.926 


0.237 


1.072 


0.259 



TABLE II. Comparison of the results for the vector mesons for the three different parameter 
sets for the effective interaction, using all eight BS amplitudes (top), and using the five leading BS 
amplitudes only (bottom). 
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luli angular 


P 




K 




J. 




calculation 


nip 


Jp 


niK* 


JK* 




u 


all 8 amplitudes 


n 7/10 
U. (4z 


U.zUr 


U.yoD 


U.z41 


1 n^o 


u.zoy 


Fi only 


U.88 


r\ on 
U.zU 


> 1.2 




1 O /I 

1.24 


n or\ 
0.20 


i^i, I'2, and i<3 


n c\r\ 

U.9U 


0.17 


> 1.2 




1 o c 


O.zO 


ri, i<A, and 


U. 1 11 


0.z3 


o.yii 


O.ZD 


1.059 


n oo 
O.ZO 


ti . .. i<5 


n 'yon 

u. / zy 


o.iyy 


o.yiy 


O.zo 


l.Ooz 


O.zoU 


leading Cliebysh 


ev decomposition 












all 8 amplitudes 


0.743 


0.211 


0.92 


0.24 


1.074 


0.262 


Fi only 


0.875 


0.20 


1.09 


0.22 


1.24 


0.22 


Fi, F2, andFs 


0.900 


0.17 


1.10 




1.25 


0.20 


Fi, Fi, andF5 


0.724 


0.23 


0.90 


0.26 


1.062 


0.28 


Fi . . . Fs 


0.730 


0.201 


0.91 


0.23 


1.065 


0.251 



TABLE III. The influence of the different covariants and of the angular dependence of the 
amplitudes on the vector meson properties with parameter set u = 0.4 GeV, D = 0.93 GeV^. For 
this table, we have calculated the loop for the decay constant up to the renormalization point 
fi = 19 GeV, since for some of the approximations considered this integral is ultraviolet divergent. 
In the case of a convergent integral, the error made by cutting off the integral at the renormalization 
point is less than 1%. 
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